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Oj, Abstract 

.^^ ' In this paper we present the necessary and sufficient conditions of separa- 

bility for multipartite pure states. These conditions are very simple, and they 
don't require Schmidt decomposition or tracing out operations. We also give a 
necessary condition for a local unitary equivalence class for a bipartite system 
in terms of the determinant of the matrix of amplitudes and explore a variance 
^ , as a measure of entanglement for multipartite pure states. 
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O ■ 1 Introduction 



Dh: 



Notation: M+ is the complex conjugate of transpose of M. 

Let \ip) and |0) be two pure states of a composite system AB possessed by 
r^ ■ both Alice and Bob, where system A (B) is called Alice's (Bob's) system. By 

C^ I Nielsen's notation \i(j) ~ \(j)) if and only if \ip) and \(j)) are locally unitarily 

^ ■ equivalent HJ. Let p4 and p4 be the states of Alice's system. It is known that 

1^) ~ |(/)) if and only if p4 and p^ have the same spectrum of eigenvalues ^ 
|2] . A pure state is separable if and only if it can be written as a tensor product 
k> , of states of different subsystems. It is also known that a state \ip) of a bipartite 

^ • system is separable if and only if it has Schmidt number 1 ^. Clearly it is 

^ ' essential to do Schmidt decomposition to find the eigenvalues of p4 and p4. To 

obtain a Schmidt decomposition of a pure state \'ip)^ we need to compute (1) 
the density operator p^^; (2) the reduced density operator p^ for system A; (3) 
the eigenvalues of p4. However it is hard to compute roots of a characteristic 
polynomial of high degree. 

Peres presented a necessary and sufficient condition for the occurrence of 
Schmidt decomposition for a tripartite pure state ^ and showed that the pos- 
itivity of the partial transpose of a density matrix is a necessary condition for 
separability 0. Thapliyal showed that a multipartite pure state is Schmidt 
decomposable if and only if the density matrices obtained by tracing out any 

^The paper was supported by NSFC (Grant No. 60433050), the fundamental research fund 
of Tsinghua university NO: JC2003043 and partially by the state key lab. of intelligence 
technology and system 



party are separable [HI ■ In [7] the local invariants of quantum-bit systems were 
investigated. In [S] [5j the local symmetry properties and local invariants of pure 
three-qubit states were discussed, respectively. In JU] the classification of three- 
qubit states was given. Bennett reported measures of multipartite pure-state 
entanglement in ^J. Meyer and Wallach J2] proposed a measure of n— qubit 
pure-state entanglement. Nielsen used the majorization of the eigenvalues of the 
reduced density operators of a composite system AB to describe the equivalence 
class under LOCC transformations. 

For a multi (n)— partite system, in this paper we illustrate the reduced den- 
sity operators obtained by tracing out the ith subsystem pi2-(»-i)(»+i)-n — 
tri^p^"^-"-) = MiM^, where i = l,2,...,n and M^ are the d"^^ x d matrices, 
of which every entry is an amplitude of the state in question. For a bipartite 
system AB, the reduced density operator p^ (p^)= AfM+, where M is the 
matrix of the amplitudes. Hence det(p4) = | det(M)p. However, for a multi 
(n)— partite system. Mi are not square. In section 2, we present a necessary and 
sufficient condition for separability for a bipartite system in terms of the deter- 
minants of all the 2x2 submatrices of the matrix of the amplitudes. Section 
3 contains three versions of the necessary and sufficient separability criterion 
for a n— qubit system. Section 4 is devoted to study the separability of multi- 
partite pure states, and two versions of the necessary and sufficient separability 
criterion are proposed. Section 5 gives a simple necessary criterion for |^) ~ |0) 
for a bipartite system. Section 6 suggests an intuitive measure of multipartite 
pure-state entanglement. 



2 The separability for a bipartite system 

Let |'0) be a pure state of a composite system AB possessed by both Alice and 
Bob. In this section we give a simple and intuitive criterion for the separability. 
Let \i) {\j)) be the orthonormal basis for system A [B). Then we can write 

IV") = Y.i,j fly l«)|j>> where YJij'=o ky P = 1- Let M = (aij)nx« be the matrix 
of the amplitudes of jV')- Then the criterion for the separability is as follows. 

1-0) is separable if and only if the determinants of all the 2x2 submatrices 
of M are zero. 

This criterion for the separability avoids Schmidt decomposition. To com- 
pute the determinants, it needs n^(n — 1)^/2 multiplication operations and 
n?{n— 1)^/4 minus operations. 

Proof. Suppose that systems A and B have the same dimension n. By 
definition, \il}) is separable if and only if we can write j-^) = (X^T^To ^^il*)) ® 
(E"=o J/jb')): where YJ-I^ |a;,p = 1 and = TJjZl |yjf = 1- By tensor product 
|'0) = ^j ■ Q Xiyj\i)\j). It means that |^) is separable if and only if Xiyj = aij, 

i,j = 0,l,...,(n- 1) (1). Let TO = ( "'' "''' ) be any 2 x 2 submatrix 

V aji ajk ) 
of M. It is easy to check det(TO) = aiia.jk - ai^aji = XiyiXjijk - XiykXjyi = 0. 
Therefore if \ip) is separable then the determinants of all the 2x2 submatrices 



of M are zero. 

Conversely, suppose that the determinants of all the 2x2 submatrices 

/ ^0 \ 

of M are zero. We can write M in the block form, M = 

\ A„_i ) 

{Bq, Bi, ..., Bn-i), where Ai is the ith row and Bi is the ith column of M, respec- 
tively, I = 0,l,...,(n- 1). Let |x,|^ = A^A+ (2) and \yj\^ = B+B, (3), 

i,j = 0, l,...,(n — 1), respectively. Under the supposition we can show that 
the above Xi in (2) and yj in (3) satisfy (1). Let us consider the case in which 
all the aij are real. It is not hard to extend the result to the case in which 
all the aij are complex. We only show |a;oj/o| = |aoo| and omit the oth- 



ers. From (2) and (3), |a;oyo|' = A^A+B+B^ = (EjJo |ao,|')(Ero' |a«l') 
Z]"7=o I'^O:;! Ifliol = I]I!7=o I "00 1 l"*jl ^ I "00 1 • In the last but one step we 
use the equality jaojl |aio| = |aoo| \o-ij\ , which holds since ■' I is 

a 2 X 2 submatrix of M. This completes the proof. 
Corollary 
If IV') is separable then det(M) — 0. 

3 The separability for a n— qubit system 

Let IV') be a pure state of a n— qubit system. Then we can write \ip) = 
J2ii,i2,....i„e{o,i}'^iii2-ir.\'ii'i2—in)- Let the density operator pi2-" = |i/')(V'l and 
^i2...(i-i)(i+i)...Ti j^g |-]-^g reduced density operator obtained by tracing out the ith 
qubit. Then pi2...(»-i)(*+i)...« = ^^.(^i2...n) ^ MiM+ , where i = l,2,...,n and 
Mi are 2^""^) x 2 matrices of the form {ab^b2...b,_i0b,+i...b„,abib2...b,_iib,+i...b„) 
in which 61,62, .-.,&« G {0, 1}. 

For example, let \ip) be a state of a 3-qubit system. Then \tp) can be written 

/ ao ai \ 



as \tp) — X]i=o"«l*)- ^3 is a 4 X 2 matrix 



Each entry of M3 is 



02 03 

\ ae ar J 
an amplitude of li/i). 

There are three versions of the separability. 

Version 1. |V') is separable if and only if the determinants of all the 2x2 
submatrices of Mi,M2,.... and M„ are zero. 

The proof of version 1 is similar to the one for a bipartite system in section 
2. 

Version 2. \ip) is separable if and only if aiUj — akai, where i + j ^ k + I 
and i (B j = k (Bl where < i,j,k,l < 2" — 1 are n— bit strings and ® indicates 
addition modulo 2. 

For example, 2, 7, 5 and 4 can be written in binary numbers as 010, 111, 101 
and 100, respectively. It is well known 010 + lll(modulo 2)= 101, 101 + 



100 = 001(modulo 2). Therefore 2+7 ^ 5+4(modulo 2) though 2+7 = 5+4 = 9. 

Using this condition it is easy to verify that states |M^) = l/-\/n(|2°)+ |2^) + 
... + |2"-i) and \GIiZ) = 1/V2(|0(")) + |l("))) for a n-qubit system [Tg are 
entangled. 

Let i\ii...im jiJ2---Jm kik2---kn and lil2-.-ln be ti— bit strings of i,j, fc and Z, 
respectively. Then version 3 is phrased below. 

Version 3. |'0) is separable if and only if aiUj — akai, where {it,jt} = {kt, h}-, 
t = l,2,...,n. 

The following lemma 1 shows that versions 2 and 3 are equivalent to each 
other. 

Lemma 1. i + j — k + I and i (B j — k Q) I ii and only if {it,jt} = {kt, h}, 
t = l,2,...,n. 

The proof of lemma 1 is put in appendix A. 

We argue version 3 next. 

Assume that jV') = (4^'|0> + 4^^|1)) ® (4^^|0) + xp^|l)) (g) ... ® (4"''|0) + 
x^ \1)). By tensor product x\ xl ....x^ — ai^i^,,^^^ where it = 0, 1, t = 

1 o rp, (1) (1) (2) (2) In) (n) , (1) (1) (2) (2) (n) (n) 

1,2, ...,n. Then a^aj = x\^' x)^' x\J x)J ....x\J x)J anAauai ^ xl>x\J xl^x\J ....xiy\J . 
Explicitly, a^aJ = auai whenever {iujt} = {kt-,lt}, i = 1,2, ...,n. 

Conversely, suppose that Qittj — a^ai whenever {it,jt} = {kt,lt}, t = 



1, 2, ..., n. Let 
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= E»i,..,»,_i,»,+i,..,»„6{o,i} |an«2,...,»J , where i = 1, 2, ...,n. 



We can show \x\_^ x\K...x]^^^ — \aiyi^...iS' ■ We only demonstrate the cases of 
n — 2 and 3 to give the essential ideas of the general case. 

When n = 2, see section 2. When n = 3, see appendix B. The two cases sug- 
gest that it be simpler to prove \xl\x^ ....xY^\^ = |aiii2...i„P (S |aiii2...i„P)" • 
Now we finish the argument for the real number case. It is not hard to extend 
the result to the complex number case. 

4 The separability for a multi (n)— partite sys- 
tem 

Assume that each subsystem has the same dimension d. Let \it] be the or- 
thonormal basis |0),|l),...,|((i — 1)) for the iih subsystem. Then any pure state 
IV') can be written as 1-0) = Z]ii7i2,...,i„=o '^*i«2---i„Ni*2.-.in)- Assume that \'\\}) 
is separable. Then we can write \i)) = {y^^'l^^i^Vi)) ® (Z]t='o ^^If N2) j ® 

■■■ ® (Si'io^lr I*"))- -^y tensor product x\^ x\^ ....x^J = aiji^...^^, where 
ii,i2,...,in e {0,l,...,(d- 1)}. 

Let the density operator p^^-" = |'0)('0I and p^^''-*^^-"^'^^^'" be the reduced 
density operator obtained by tracing out the ith subsystem. Then pi2-(«-i)(2+i) ■■" — 
tri{p^'^--'^) = MiM^, where i — 1,2, ...,n and Mi are d"^^ x d matrices of the 
amplitudes of the form 

(afcife2...fc.-i0fei+i...fc„,afcifc2...fci_iife,+i. ..*;„, ■•■,afeife2...fc,_i(d-i)fc,+i...fc„), where A:i,fc2,...,fcj_i,/i;i+i,...,/c„ e 
{0,l,...,(d-l)}. 



There are two versions of the separabihty. 

Version 1. |V') is separable if and only if the determinants of all the 2x2 
submatrices of Mi, M2, ... and M„ are zero. 

Version 2. Jt/)) is separable if and only if aiji2...i„ajjj2...j„ = akik2...k„ai^i^,,,i^, 
where {it,jt} = {hJt}, t = l,2,...,n. 

The proof of version 1 is similar to the one for a bipartite system. The proof 
of version 2 is similar to the one for a n— qubit system. 

When n = 2, the criterion is reduced to the one for a bipartite system. When 
d = 2, the criterion is reduced to the one for a n-qubit system. 

5 A necessary condition for a local unitary equiv- 
alence class for a bipartite system 

We use the following lemma 2 to establish the necessary condition. 

Lemma 2. Let \tp) be a pure state of a composite system AB possessed by 
both Alice and Bob. Let M = {ajk)nxn be the matrix of the amplitudes of IV')- 
Let p^^ = \ip){ip\ and p"^ = trsip^^)- Then |det(M)|^ is just the product of 
the eigenvalues of p^. 

The proof is put in appendix C. 

Lemma 2 reveals the relation between the determinant of the matrix of the 
amplitudes and the eigenvalues of p"^ for a bipartite system. 

The corollary of lemma 2 

Let M^ {Mff,) be the matrix of the amplitudes of a pure state \il}) {\4>)) of a 
composite system AB. Then |det(M^)| = |det(M0)| whenever jV') ~ \4>)- That 
is, |det(M^)| is invariant under local unitary operators. 

It is well known that it only needs 0{tt') multiplication operations to com- 
pute |det(M)| instead of doing Schmidt decomposition in P|5]. 

For a two-qubit system, let \ip) = a|00)+6|01)+c|10)+d|ll) and p^^ ^ |-0)(^|. 

By lemma 2 \ad — 6cp is the product of the eigenvalues of p^ . Let \ad — hc\ —G. 
We can show that e satisfies <G< ^ and eigenvalues \± = ii^^-^=^^. Hence, 

6 The variance as a measure of entanglement 

We obtain the necessary and sufficient conditions of separability in sections 
2, 3 and 4. Apparently, \at^t2...i„aj^j^,,,j^ - akik2...k„ai^i2...ij, where {itjt} = 
{kt, It}, i = 1, 2, ..., n., is just a deviation from a product state. It is intuitive to 
suggest the variance: J2 \0'ili2■■■i„ajlj2...j,^ - akik2...k„ai^i2...ij , where {it,jt} = 
{ktjt}: i = 1,2, ...,n, as a measure of entanglement of \ip). Let De{\'(P)) be the 
measure of entanglement. 

De{\^)) has the following properties. 

Property 1. DE{\tp)) = if and only if jV') is separable. 

The properties for a two-qubit system 



For a two-qubit system, let |V') = a|00) + 6|01) + c|10) + rf|ll). Then 



DE{m^\ad-bc\ 



2 



Property 2. The maxhuum of DeHiP)) = \ad-bc\ < {\ad\ + |6c|)^ < 



1 



V 2 "•" 2 -' 4- 

When a, 6, c and d are real, by computing extremum it is derived that the 
maximally entangled states must be of the forms: x|00) + y|01) — y\10) + 2;|11) 
or a;|00) + y|01) + y|10) - a;|ll). 

Property 3. |V') -- |V^') if and only if DEQi^)) = DeHtP')). 

Given \ij}) = a|00) +6|01) +c|10) + d|ll) and \ip') = a'|00) +6'|01) +c'|10) + 
d'|ll). Suppose that lip) ~ \ip'). By the necessary condition in section 5, 

DEim ^ DEim). 

Conversely, suppose De{\^)) = De{\iP'))- Let us show lip) ^ lip'). Using 
Schmidt decomposition, we can write \ip) ^ -v/Ai|00) + V^IH); where Ai + 
A2 = 1. As discussed above \ad — bc\ = VAiV^- As well using Schmidt 
decomposition we can write ji/i') ~ ^/pTlOO) + ^/p2\ll), where pi + P2 = 1, and 
\a'd'-b'c'\^^^. ThusAiA2-piP2. Then Ai(l - Ai) ^ pi(l-pi). There 
are two cases. Case 1. Ai = pi. Then A2 = P2- Case 2. Ai + pi + 1 = 0. In the 
case A2 = pi and Ai — p2- It means that \ip) and \'ip') have the same Schmidt 
co-efRcient for either of the two cases. By factor 5 in PPI|2], \ip) ~ |'0')- 

Nielsen in ^J showed |^') ~ \'4>") by calculating eigenvalue, where ji/'') = 
^/^|00) + 7ar|ll),and|V'") = (|00) + |l)(cos7|0)+sin7|l)))/\/^. By property 
3 it only needs to check JctJ^Jot- = sin 7/2. 

Conclusion 

In this paper we have presented the necessary and sufficient conditions of 
separability for multipartite pure states. These conditions don't require Schmidt 
decomposition or tracing out operations. By using the conditions it is easy to 
check whether or not a multipartite pure state is entangled. 

Appendix A. The proof of lemma 1 

Let aia2---an, f3iP2---Pn, Sid2---Sn and 71 72-. .771 be the n— bit strings of a, 
P, 6 and 7, respectively. 

Lemma 1. {a^, (3i}—{di, 7^}, i — 1,2, ..., n, if and only ifa + /3 = (5 + 7 and 
a(B /3 = S (B 7, where ® indicates addition modulo 2. 

Proof. Suppose {at. Pi}— {Si, ji}, i — l,2,...,n. Since a + f3 = (ai + 
/3i)2"-i + (a2 + /32)2"-2 + ... + (a„ + /?„) and <5 + 7 = (<5i + 7i)2"-i ^ (^^ ^ 
72)2""^ + ... + {Sn + 7n), by the supposition it is easy to see a + P — S + j. It 
is straightforward to obtain aia2...a„ © PiP2---Pn = SiS2---5„ ® 7i72...7n- 

Conversely, suppose a + P — S + j and a® P — S (B j. First let us consider 
the case where n — 1. There are three cases. 

Case 1. ai + Pi — Si + ji = 0. This means ai ~ Pi ~ Si ~ 'fi ~ 0. 

Case 2. ai + Pi = Si + ji = 1. This implies {ai,Pi} = {^1,71} — {1,0}. 

Case 3. ai + Pi = Si + ■^i = 2. This says ai = Pi = Si = ji = 1. 

No matter which of the above three cases happens, it yields {ai. Pi}— {Si, 

Til- 
Let us consider the case n. Since a + P = S + ■y, (ai + /3i)2"^^ + (a2 + 

/32)2"-2 + ... + (a„+/3„) = (5i+7i)2"-i + ((52+72)2"^2 + ... + (^„+7«)- Again 



since a ® /? = (5 ® 7, that is, aia2...an ® /3i/32---/3n — SiS2---Sn 7i72---7n, we 
obtain ai (B Pi = Si (B 7i, i — 1,2, ..., n. There are two cases. 

Case 1. a„ e /3„ = 5„ e 7„ = 1. In the case {a„, /3„} = {(5„, 7„} = {0, 1}. 
Then (ai + f3i)2"'^ + (^2 + /32)2"-3 + ... + (a„-i + /3„-i) = ((5i + 71)2"-^ + 
(52+72)2"~^ + ... + ((5„-i+7n-i) and a,® A =(5,®7„ i = l,2,...,n-l. By 
induction hypothesis {ai,Pi} — {Si,ji}, i — l,2,...,n— 1. 

Case 2. «« ® /3n = <^ti ® 7n = 0. There arc two subcases. 

Subcase 2.1. a„ = /3„ = (5„ = 7„ = or «„ == /3„ = (5„ = 7„ = 1. As 
discussed in case 1, we can obtain {ai,Pi} — {Si,"fi}, i = l,2,...,n — 1 by 
induction hypothesis. 

Subcase 2.2. a„ = /3„ = 1 and (5„ = 7„ = or a„ = /3„ = and Sn — Jn — 1- 
Let us consider the former case. In the case (ai + /3i)2"^^ + (a2 + /32)2"^^ + 
... + (a„_2 + /3«-2)2 + (a„-i + /3„-i + 1) = 

{Si + 7i)2"-2 ^ (^2 ^ 72)2"-3 ^ _ ^ (^^^_2 + 7«-2)2 + (<5„_i + 7„_i). 

Since a„_i ® /3„-i = (5ri-i ® 7ri-i, either a„-i ® /3„-i = (5„_i ® 7„_i = or 
1 causes that one of (a„_i + /?„_i + 1) and {Sn-i + 7n-i) is odd and the other 
is even. It contradicts a ® /3 = (5 ® 7. 

Appendix B. The separabihty for a n— qubit system 

When n — 3, let us show |a;^^ x^-^ x^-^ ^ = loijijigp when aiaj — auai, where 
{it, jt} = {kt,lt}, t = 1,2,3. We only illustrate |4^^x[,^^x^^^|2 ^ |^^^^|2_ Q^j^g^. 
cases then follow readily. Experientially, it is simpler to prove |xq Xq Xq p = 
|aoooP(Ejj,fce{o,i} l««jfcH(E»,j,fce{04} |a»jfcP), where |x[,^^|2 = E»jg{o,i} I^o^I^' 



X, 



(2)12 _Y- U,„,|2 ,„j u(3)|2 _ Y- u „|2 







Efc,iG{o,i} '"'^o'l^ ^'^'^ 1^0 P = Ep,ge{o,i} I"P«"I 



First we show that aoijakmcipqo can be rewritten as aoaoO'aia2a3(^Si52S3- There 
are the following four cases. 

Case 1. Consider aoijUkoi and the pairs {0, fc}, {i, 0} and {j, I}. If j * Z = , 
then aojjafeoz = aoooOfciO+o since {j, ?} = {0, j + ?}. 

Case 2. Consider aoijUpgo and the pairs {0,p}, {i, g} and {j, 0}. If i * q — 0, 
then aoijUpqa = aoooap(i+g)j since {i,g} = {0,i + q}. 

Case 3. Consider akoia-pqo and the pairs {fc,p}, {0, q} and {/, 0}. If fc *p = 0, 
then akoiapqo = aoooa(k+p)qi since {k,p} = {0, fc + p}. 

Case 4. Otherwise i=j~l~k^p = q = l. It is not hard to derive 
030506 == 010705 — OoO^. 

Second, let us show that 00000^102030515253 can be rewritten as aoijakoiCLpqo- 
If OoooOQiQ2a305i52<53 is of the forms: ooooOoyOfeo;, ooooOoijOpgo or ooooOfeo;OpgO) 
then these forms are desired. Otherwise ooooOaiQ2Q30(5i52i53 must be aoOgag, 
O0O3O3, 000505 or of the form 00070^^4, which can be rewritten as 020405, O1O2O3, 
010405, oiOoOrst, respectively. 02O4O6, 01O2O3 and 01O4O5 are just desired and 
oiOgOrst is furthermore rewritten as follows. There are three cases. 

Case 1. In the case r = or s = 0, this is desired. 

Case 2. In the case r = s = t = 1, 010307 = 030500, desired. 

Case 3. In the case r = s = I and t = 0, oiogog = 020500, desired. 

Appendix C. The proof of lemma 2 

Proof. Suppose that systems A and B have the same dimensions n. Let \^p) = 



YJij=oO'i3\i)\i)- ThcnM=: (ay)„x«- Let density operator p-^-^ = |i/')(V'l- Then 

p-'^ - (Epo«d*>b-»(E".~io«r.(^i(^i) = Y:~l,Y.rki,a,,a'i,\m{im 

= Eri=o E?fe=o '^y'^rfcK)b)(^l(^l- The reduced density operator for system 
A is defined by p^ — trsip^^)- Let us compute p^. 

P'^ = Y,1J^oY,'^jj^^Qaija*ik\i){A^kj (where Stj = 1 when k ^ j. Other- 
wise 0.) = E"rioEj=o ay-ajj-NX^I = E"rio(Ej=o aya?j)N>(^|- Let A, = 
{aio,aii, ....ai(„_i)), that is, the ith row of A. Then E?=o ^ij^tj ~ ^^i^i ■ Fi- 



nallyp^^E-rio^^A+IOai 



/ ^0 \ 
Ai 



{A+,A+,...,A+_^) ^ MM+. Thus 



det(p'^) — I det(M)p. Hence |det(M)| is just the product of the eigenvalues of 
p^. Q.E.D. 
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